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1 INTRODUCTION 



ABSTRACT 

We examine the alignment between Brightest Cluster Galaxies (BCGs) and their host 
clusters in a sample of 7031 clusters with 0.08 < z < 0.44 found using a matched-filter 
algorithm and an independent sample of 5744 clusters with 0.1 < z < 0.3 selected 
with the maxBCG algorithm, both extracted from the Sloan Digital Sky Survey Data 
Release 6 imaging data. We confirm that BCGs are preferentially aligned with the 
cluster's major axis; clusters with dominant BCGs (> 0.65 mag brighter than the mean 
of the second and third ranked galaxies) show stronger alignment than do clusters 
with less dominant BCGs at the 4.4cr level. Rich clusters show a stronger alignment 
than do poor clusters at the 2.3a level. Low redshift clusters {z < 0.26) show more 
alignment than do high redshift {z > 0.26) clusters, with a difference significant at 
the 3. Oct level. Our results do not depend on the algorithm used to select the cluster 
sample, suggesting that they are not biased by systematics of either algorithm. The 
correlation between BCG dominance and cluster alignment may be a consequence 
of the hierarchical merging process which forms the cluster. The observed redshift 
evolution may follow from secondary infall at late redshifts. 

with their host clusters from SDSS data in this paper. 

BCGs are not simply the brightest galaxy picked from 



The present cosmological paradigm ( Komatsu et al. 2009 1 



the luminosity function of a cluster. Humason et al. ( 1956 1 



predicts the existence of coherent structures in the universe 
as well as alignments between structures on various scales. In 
a cold dark-matter dominated universe with structures form- 
ing hierarchically from the bottom up, filamentary struc- 
ture formation and tidal forces may cause alignments be- 



Tremaine & Richstone (19771; Lauer (19881; Postman & 



Lauer (19951; Garijo et al. (19971; Dubinski (19981; Loh & 



Strauss (20061 and others point out that BCGs have prop- 



tween different elements of the hierarchy (e.g. Torlina et al. 
|2007[ [Faltenbacher et al.|2008[ ). Simulations using a ACDM 
cosmology have predicted alignments between clusters and 



super-clusters (e.g. Basilakos et al. 20061, alignments be- 



tween clusters ( [Hopkins et al.|2005[), an d alignments between 



galaxies and clusters (Dubinski 19981. Observational stud- 



ies have found alignments between galaxies and large-scale 
structures (IFalte nbacher et al.||2009t ; clusters and clusters 



(e.g..|Binggcli 19821, galaxies and clusters (e.g.. 



BinggeU||1982j |Struble||1990l |Kim et al.|[2002 



Sastry 



1968 



erties that suggest a formation scenario distinct from that of 
other galaxies, including photometric and colour homogene- 
ity, a metric luminosity independent of cluster richness and 
significantly higher than that of lower-ranked galaxies, dis- 
turbed morphologies, extended low surface brightness stellar 
haloes and central location. 

The shapes of clusters, the special properties of bright- 
est cluster galaxies, and their alignment are clues to the 
process by which clusters formed. Understanding the align- 
ment offers an observational test of the current cosmological 
paradigm. The alignment of BCGs with their hosts was first 



Hashimoto 



et al. 12008 1, and galaxies within groups (e .g.. |Yang et al 
2006[ [Fahenbacher et aL]|2007[ [Wang et al"1|2008[ ). Most of 
the earlier observational probes of galaxy-cluster alignment 
were based on small samples and only have shallow redshift 
coverage. With the availability of large imaging surveys such 
as the Sloan Digital Sky Survey (SPSS , York et al.t20 00) we 
can now extend these earlier observations. In this spirit we 
study the alignment of brightest cluster galaxies (BCGs) 



(19801. Detailed observations were carried out by Binggeli 



(1982 



noted by Sastry ( 1968 1 and studied by Carter & Metcalfe 



(hence the alignment is also known as the Binggeli 



effect); Struble & Peebles ( |1985 1 (who found no alignment 



with a small cluster sample 
govski (|1987[ ) ;|Lambas et al 



et al. 



et al. 



( 1992 1; Fuller et al. ( 1999 1; Kim et al 



Rhee fc K atgcTF ("1987") ; |Djor 
1988|); [Struble (,1990^ ; |Trevese 



(20061 and Siverd et al. (20091. The effect is observed 



(2002l;ponoso 



when the cluster shape is defined by the distribution of mem- 
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ber galaxies or by the observed distribution of x-ray emitting 
hot gas ( Hashimoto et al.|2008 l. The later studies were able 
to confirm the alignment for clusters with redshifts up to 
z = 0.5. The effect is strongest for cluster shapes measured 
for red, centrally concentrated galaxies, and it weakens with 
increasing separation between the central galaxy and the 
satellites. 

The two leading candidate mechanisms to explain the 
alignment of the BCG with their parent cluster are that both 
are formed from infall along preferred directions in filaments 
(e.g. Dubinski|1998 1, and that BCGs are aligned by tidal in- 
teractions (e.g. Faltenbacher et al.|2008 l. The points in cos- 
mic history at which these mechanisms act are very different: 
whereas filamentary infall is likely to affect alignments dur- 
ing and immediately after cluster virialization tidal affects 
can act during the cluster's entire lifetime. |Ciotti fc Dutta| 
1 1994 1 have shown that the time scale on which a prolate 



galaxy's orientation is affected by a cluster's tidal field is 
much shorter than a Ifubble time. Present-day alignments 
may therefore either be the result of primordial alignments 
stemming from the period of cluster formation or can have 
grown during the cluster's lifetime. By studying the redshift 
evolution of the alignment effect we can hope to distinguish 
these two cases. 

In this paper we use Sloan Digital Sky Survey Data 
Release 6 data ( York et al.|200d Adelman-McCarthy et al 



2008) to study the alignment effect in 12,755 clusters ex- 



tending out to z = 0.44. In §2 we describe our cluster selec- 
tion and in §3 we analyse the dependence of alignment on 
BCG dominance, cluster richness, and redshift. We discuss 
the implications of our findings in §4 and summarise in §5. 
Throughout this paper we assume an Qm = 0.3, JIa ~ 0.7, 
and Ho — 70 km s~^ Mpc~^ cosmology. 



2 DATA 



2.1 The Sloan Digital Sky Survey 



The Sloan Digital Sky Survey (SDSS; [York et al.|[2000l ) is 
an imaging and spectroscopic survey that has covered one- 
quarter of the celestial sphere. It is carried out with a dedi- 
cated 2.5m tel escope (|Gunn et al.|2006| ) using a large format 
CCD camera ( Gunn et al.|1998 |. The imaging data are col- 
lected in drift-scan mode in five band passes {u,g,r,i,z) with 
effective wavelengths of 3551, 4686, 6165, 7481, and 8931 
A. The imaging data are automatically reduced through a 
series of software pipelines which find and measure objects 



and provide photometric and astrometric calibrations ( Lup- 
ton et al.|2007| [Lupton et al. 2002 ; Pier et al.||2003| [Tucker 



et al. 20061. The photometric calibrations are accurate to 
about 1% rms in g, r, and i, 3% in u and 2% in z for bright 
(< 20 mag) point sources (Ivezic et al. 20041. We restrict 



ourselves to the objects that have reliable photometric data 
by using the SDSS clean photometry flags. Targets for spec- 
troscopy are selected from the imaging data on the basis of 
their photometric properties. A pair of dual fibre-fed spec- 
trographs ( [Uomoto et al.|[l999[ ) can observe 640 spectra at 
a time with a wavelength coverage of 3800 A to 9200 A and 
a resolution of 1800 A to 2100 A. 

The SDSS data are described in the data release 
papers (Adelman-McCarthy et al. 2008 for the Sixth 



Data Release which we use) and are documented at 
http:/ /www. sdss.org. 

2.2 Cluster Catalogues 

In this paper we use two catalogues of clusters of galaxies. 
One by Koester et al. ( 2007a|b I takes advantage of the con- 
centration of cluster galaxies in colour-magnitude space on 



the red sequence, and the other, by Dong et al. (20081 iden- 



tifies clusters by matching galaxy distributions in position- 
magnitude space to an a priori filter. The two catalogues are 
likely to have different systematic errors such as miscenter- 
ing and problems caused by cluster overlap, so computing 
results from the two allows us to reduce the severity such 
problems are. 

There is a well known bimodality of galaxies in colour- 



colour space (e.g. Strateva et al. 2001 Baldry et al. 2004 



Bell et al. 2004, Cassata et al. 



2008) with red quiescent 
spheroidal galaxies occupying what is known as the red se- 
quence, and blue actively star-forming galaxies concentrated 



in the blue cloud. As first noted by Baum (19591, red se- 



quence galaxies are predominantly found in clusters of galax- 
ies. Since the location of the red sequence in colour-colour 
space is dependent on redshift via the K-correction it is a 
viable photometric redshift indicator for galaxy clusters and 
field galaxies (e .g..'Anni s et al.|1999}|Gladders fc Yee|2000 



[Koester et al.| |2007a b; |Eisenhardt et al.[[2008[ and refer- 
ences therein). These properties make the red sequence an 
ideal tool for identifying galaxy clusters in surveys (e. g..|Yee 
eraL][T999l [Gladders fc Y ce 200 0] [Lubin et al.[[205ol [Glad 



ders fc Yee[[2005 Koester et al.| 2007a I . The red sequence 



allows one to isolate galaxies at a given redshift and thus 
remove contamination by projected foreground and back- 
ground galaxies. 

[Koester et al.| (|2007a|b| ) follow [Gladders fc Yee| ( [2000 \ 
and[Annis et al. ( 1999 1 using a technique known as maxBCG 



to search directly for over-densities of red sequence galax- 
ies. They take advantage of the uniform colours and lumi- 
nosities of BCGs in the selection procedure. The likelihood 
that each galaxy in a photometric sample has the photo- 
metric properties of a BCG and resides in an over-density of 
red sequence galaxies is evaluated at a grid of assumed red- 
shifts. The redshift which maximises the likelihood (hence 
the name maxBCG) is used as a first estimate of the cluster 
redshift. Using this centre, all red sequence galaxies within 
1 h~^ Mpc projected radius are potential cluster members. 
A percolation technique is then used in order to determine if 
the galaxy in question is in fact the cluster centre. The cat- 
alogue sample we employed contains 12, 766 clusters with 
photometric redshifts between 0.1 and 0.3 and is approxi- 
mately volume limited. Koester at al. use the maxBCG al- 
gorithm on realistic mock catalogues and find that the sam- 
ple is more than 90% pure and more than 85% complete 
for clusters with masses ^ 1 x 10^* Mq. We extract BCG 
positions (which are defined to be at the cluster centres in 
this catalogue) and redshifts in order to select red sequence 
galaxies from the SDSS. Koester et al. find an rms scatter 
of only a = 0.015 between the clusters' photometric red- 
shift estimates and spectroscopic redshifts of the BCG for 
the 7813 clusters in their catalogue with spectroscopy in the 

SDSS^ 

Dong et al. ( 2008 1 use a variant of the matched fil- 
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ter technique previously used by Postman et al. (19961; 
Kawasaki et al. ( 1998 ) ; Schuecker & Boehringer ( 1998 1 ; Kep- 



ner et aL| (|1999j); ,Bramel et al.| ( |2000| ; |Kim et al.| ( |2002[ ), 
among others, which fits the distribution of galaxies in mag- 
nitude and position space to a standard spatial profile (such 
as that predicted by Navarro et al.|19^ I , using prior knowl- 
edge of spectroscopic or photometric galaxy redshifts and 
the gala:xy luminosity function. The technique does not ex- 
plicitly fit for the red sequence to select galaxies and hence 
can identify clusters of blue galaxies (if they exist). The 
algorithm generates a cluster likelihood map in position- 
redshift space whose peaks correspond to positions where 
the matches between the survey data and the cluster filter 
are optimised. Using a realistic mock SDSS catalogue Dong 
et al. show that the catalogue is « 85% complete and over 
90% pure for clusters with masses above 1.0 x 10"/i"^Mo 
and redshifts up to z = 0.45. The errors for estimated cluster 
redshifts are typically less than 0.01 (comparable to the pho- 
tometric redshift errors of individual red-sequence galaxies) 
. The positions of clusters in the catalogue are the geomet- 
ric centres of the galaxy distributions and do not necessarily 
correspond to a single galaxy position. Dong et al. do not 
identify BCGs in their catalogue. The matched filter has a 
faint end absolute magnitude cut off at 0.4L* and hence the 
catalogue is approximately volume limited. The value of Li, 
used assumes passive evolution, which may not be valid in 
detail. 



2.3 Cluster Galaxy Selection from SDSS 

In this study we use both catalogues but will concentrate on 
the results from the Dong et al. matched filter catalogue due 
to its larger size and deeper redshift coverage. In the follow- 
ing discussion, unless stated otherwise the values quoted and 
figures shown refer to the Dong et al. catalogue. In order to 
study the alignment of the BCG with the cluster we require 
the cluster centre and redshift in order to identify member 
galaxies. We define the cluster's shape using the identified 
member galaxies and identify the BCG as the brightest of 
these members. 

To construct our sample of cluster galaxies we select 
galaxies from the SDSS database within 1 Mpc projected of 
a given search centre (the geometric cluster centre quoted 
in the catalogue). We select only red sequence galaxies at 
the cluster redshift provided in the catalogue, using red se- 
quence colour cuts g — r and r — i determined by |Loh fc| 



Strauss (20061. The model magnitudes that we use are cor- 



rected for Galactic extinction using the maps by ISchlegel] 
|et al.| ( |1998[ ). We initially performed two colour cuts on 
the data corresponding to selecting objects at a given red- 
shift that fall within the 2a and 3(j ellipses in colour-colour 
space of the Loh & Strauss red sequence. However, we found 
that our results do not strongly depend on the colour cut 
used, and hence we shall present only the data for the So- 
cut in this paper. The colour-colour ellipses used become 
larger at higher redshift to account for increasing photo- 
metric errors. One might suspect that this will also increase 
contamination by foreground and background galaxies and 
thus make the resultant clusters rounder. However, we find 
that our high redshift clusters are not systematically rounder 
than our low redshift clusters, suggesting that any redshift- 
dependent contamination is not severe. In addition to be- 



ing more cleanly separated from the field galaxies in colour 
space, and therefore showing higher contrast and less con- 
tamination, the red sequence galaxies tend to be deeper in 
the potential well of their clusters, and therefore will do a 
better job of refiecting the true cluster shape. 

After performing the colour cuts, the brightest galaxy 
in the i-band within the 1 Mpc radius from the search centre 
is selected as the BCG. In 85% of the clusters in common 
between the Dong et al. and Koester et al. catalogues, our 
identification of the BCG matched that of Koester et al. 
As BCGs are observed to be a highly uniform population 
in colour and luminosity (e.g.. Postman fc Lauer|1995 l and 
since the SDSS probes a uniform range of magnitudes at 
each redshift, we impose a cut of i < iscG — 3 on the mem- 
ber red sequence galaxies to ensure that we sample each 
cluster approximately to equal depths. At the highest red- 
shift in the Dong et al. catalogue {z = 0.44) we sample the 
cluster down to i = 21.67, where the SDSS is still quite com- 
plete ( [Yasuda et aL]|2001 l. 

The preliminary centre of each cluster is defined as the 
mean position of all selected red sequence galaxies. The final 
cluster galaxy sample consists of all red sequence galaxies 
within 0.5 Mpc of this preliminary centre and the final geo- 
metric cluster centre is taken as the mean position of these 
galaxies. We find that the ellipticity of our clusters does not 
strongly depend on using a 1 Mpc or 0.5 Mpc radius to select 
cluster members, nor are rich and poor clusters affected dif- 
ferently. This is a result of determining cluster shapes using 
radius-weighted second moments, which we describe below. 

Of the 36, 785 clusters in the Dong et al. catalogue we 
can extract data for 23, 106 clusters (the others being outside 



the redshift range for which Loh & Strauss ( 2006 1 empiri- 



cally define the red sequence). We discard clusters in which 
the geometric centre is farther than 0.5 Mpc from the origi- 
nal search centre (this affects 1395 clusters) as well as those 
clusters in which the BCG is farther than 0.5 Mpc from 
the geometric centre (7658 of the remaining clusters). In the 
maxBCG routine of Koester et al. this BCG-miscentering 
may result from the algorithm itself. In a cluster with no 
clearly dominant BCG, the choice of BCG by the cluster- 
finder is less well-defined. Since this routine defines the BCG 
as the cluster centre, by removing clusters with clearly off- 
centre BCGs we may in fact be removing clusters with less 
dominant BCGs. However, since our alignment results from 
the Koester et al. catalogue match those of the Dong et al. 
catalogue (which does not depend on the BCG to determine 
cluster centres) it does not seem likely that such a system- 
atic effect would bias our sample towards clusters with more 
dominant BCGs. 

A candidate BCG that is far from the centre may in 
fact simply be a bright foreground galaxy not physically as- 
sociated with the cluster, or may indicate that the system is 
a pair of merging clusters. In the latter case there may exist 
two BCGs that may have refiected the shape of their origi- 
nal parent cluster but not necessarily of the merged cluster. 
Figure [l] shows the histogram of projected distances of the 
BCG from the cluster centre. The BCG is within 350 kpc 
of the centre in 80% of the clusters in the Dong et al. sam- 
ple and in 90% of the clusters in the Koester et al. sample. 
Similarly [Postman fc Lauer| ( |1995l ) find that for 90% of the 
clusters in their z < 0.05 sample the BCGs lie within 350 



kpc. Skibba et al. (20101, using a galaxy group catalogue 
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Figure 1. Histogram of projected distances of the BCG from the 
cluster centre. The Dong et al. sample is shown in the solid line 
and the Koester et al. sample in the dashed line. In the Dong et 
al. sample we find that 80% of galaxies are within 350 kpc and 
50% are within 200 kpc. In the Koester et al. sample we find 90% 
within 350 kpc and 63% within 200 kpc. 



constructed from the SDSS, find that the brightest galaxy 
does not reside at the halo centre (i.e., the position of low- 
est specific potential energy) in 25% — 40% of their sample. 
The cuts imposed on the data for both the Dong et al. and 
Koester et al. catalogues are described in Table [T] 

To compare the two catalogues we have matched cluster 
centres, requiring that their projected distance be less than 2 
Mpc and their redshifts agree to within 0.02. The sample of 
clusters we compare is restricted to the redshift range of the 
Koester catalogue (0.1 < 2 < 0.3) which corresponds to 3189 
clusters from the Dong et al. sample and 6625 clusters from 
the Koester et al. sample. Our matching routine returns a 
total of 1000 matching clusters between the two samples, 
heavily weighted to the richer systems. The geometric cen- 
tres agree (separation less than 300kpc) in 990 cases and 
the selected BCGs match in 984 cases. For clusters that do 
not have matching BCGs, the initial search centres are on 
average more than 300 kpc apart. The richness distributions 
are very similar in both catalogues. 



3 ANALYSIS 

3.1 Cluster Shapes 



We follow the method used by Kim et al. ( 2002 I to calculate 



the orientation and elongation of clusters. We first calculate 
the radius-weighted second moments (the first moments 
vanish by definition of the cluster centre) for all cluster 
members: 



xy\ 
r2 / 



(1) 



where x and y are the respective distances of a given 
member galaxy from the cluster centre defined in the 
equatorial coordinate system. Using the relation between 
the Stokes parameters Q and U, the position angle {(j>) and 
axis ratio (a = ^, where a and b are the semi-major and 
semi-minor axis of the ellipse) we may determine </!> and 
a for each cluster. In order to compare the position angle 
of the clusters to those of galaxies reported in the SDSS 
database we convert 6 from West of North to East of North. 



1 — a 

Q = cos(20) = M^^ - Myy = 2M^^ - 1 (2) 

1 -|- Q 



U ; 



1 - a 
H- a 



sin(2<ji) = 2M^ 



(3) 



The uncertainties in Q and U are given by Poisson statistics 
as: 



N{N -1) 



iv(iv-i) 



(4) 
(5) 



From the above it is straightforward to solve for and a: 



D 



D 



[/2 



(6) 



The position angle and ellipticity of the BCG are provided 
in the SDSS data, based on a two-dimensional fit to a 



PSF-convolved de Vaucouleurs profile (Stoughton et al 



|2002] ). 

The cluster and BCG position angles are distributed 
isotropically, as is expected. Clusters are elliptical rather 
than round ( Figure [2| which is i n agreement with observa- 



tions by e.g.. Carter & Metcalfe ( 1980fc and Binggeh (19821 



and supported by simulations (e.g.. [Splinter et al 



1997 



[Hopkins et al.||2005| ). Our population of clusters is more 
elliptical than those of the Hopkins et al. simulations (who 
find a mean ellipticity of b/a « 0.7 a.t z = 0), with the 
richer clusters having a mean ellipticity closer to, but still 
more elliptical than, the values predicted by the simulation 
(Figure |2|. The shapes determined by Hopkins et al. are 
calculated using the entire cluster (all particles within the 
friends-of-friends radius, which is greater than 0.5 Mpc in 
all cases), whereas we employ only the central 0.5 Mpc. 
Hopkins et al. note that the measured ellipticities decrease 



with decreasing radius, however, AUgood et al. (20061 have 



noted that the central parts of dark matter halos are more 
aspherical than the overall cluster. It is possible that this 
effect is the cause of the difference in the values of ellipticity 
determined by us and those measured in the Hopkins et 
al. simulations. There is no correlation between the cluster 
and BCG ellipticities (Figure [sj, that is highly elongated 
BCGs show no tendency to live in highly elongated clusters. 
[Hashimoto et ah] ( |2008[ ) found the same result measuring 
the ellipticities of x-ray emission of clusters. 

Clusters with offset BCGs are potentially two merging 
clusters and the shape of the selected BCG may be related 
to the shape of its original host cluster,rather than the 
galaxy distribution of two merging clusters. Hence we might 
expect a weaker alignment between these clusters and their 
BCGs. We divide the sample roughly in two by splitting 
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Table 1. Cluster Samples 





Dong et al. 


Koestor et al. 


Sky Coverage 


6500deg2 


ft! 7500deg2 


Total Input Catalogue 


36785 


12766 


0.08 < z < 0.44 


23106 


12766 


Centrec eometric < 0.5 Mpc from Centre7„p„t 


21711 


12202 


BCG < 0.5 Mpc from CentreQeometric 


14053 


10754 


> 2.3 


9115 


7071 


-^Members ^ ^ 


8081 


6626 


b/agCG < 0.9 


7031 


5744 


Centre (BCG ^5 0.2 Mpc from CentreQ^g^^tj-i^) 


3618 


3647 


Offset (BCG > 0.2 Mpc from Centregeometric) 


3413 


2097 



at a BCG-centre distance of 0.2 Mpc and we call these the 
Centre and Offset samples respectively. Figure [4] shows the 
distribution of the number of red sequence galaxies for the 
final cluster sample. The Offset sample has slightly poorer 
clusters on average than does the Centre sample. 

Clusters with very small axis ratios (i.e. very elongated) 
are somewhat unusual objects. Sixty percent of the 212 
clusters in our final sample with b/a < 0.1 have offset BCGs 
(vs. 50% in the total sample), and 94% have fewer than 
8 members. These objects may be less virialized, and are 
certainly less rich systems. We find that excluding them 
does not qualitatively change our statistics or conclusions. 

We remove clusters that are consistent with being 
round (4938 clusters), as we will be unable to measure an 
alignment signal for such objects. Following |Kim et al.| 
( 2002 1 we define this condition as: 



+ 



sC 2.3 



(7) 



which marks 68% in a distribution with two degrees of 
freedom i.e. we require clusters to be farther than Icr away 
from D = (see Equation [6| which defines a round cluster. 
We also remove the 1034 remaining clusters that have four 
or fewer members on the red sequence, as a position angle 
defined from them will be dominated by shot noise. We re- 
move from our final catalogue those clusters whose BCGs 
are round, which we define as h/a > 0.9 since their position 
angles will be meaningless. This eliminates an additional 
1050 clusters. The final cluster catalogues are summarised 
in second part of Table [l] Throughout the subsequent anal- 
ysis we use the cluster sample labelled as hjascG < 0.9 in 
Table [l] as the Total sample and split this into the Centre 
and Offset samples. 




Figure 2. Distribution of axis ratios hja of clusters and BCGs. 
The solid line shows the ellipticities of clusters, as determined by 
our algorithm. The dash-dotted line shows the ellipticities of the 
259 clusters with more than 20 members. In these clusters the 
shape should be determined with higher signal-to-noise than in 
the poorer clusters. They are less elliptical than the total sam- 
ple of clusters. The difference results from the fact that sparse 
sampling in the poor clusters has a tendency to make them more 
elliptical (see Appendix A in AUgood et al.l2006) than more highly 
sampled clusters. The dashed line shows the distribution of axis 
ratios for BCGs given by the SDSS database. In our analysis we 
remove BCGs with b/a> 0.9 for which the position angle deter- 
mination is meaningless. 



3.2 Alignment 

With our measurements of the position angles of clusters in 
hand, we investigated a number of alignment signals between 
BCGs and their host clusters. We define the alignment of the 
BCG with the cluster as 



/^<f) = \<j)BCG — 4>ClustEr\ (8) 

and we say that the BCG and cluster are aligned for 



A</) ^ 30°, following Binggeli (19821 



We quantify the strength of the alignment signal by cal- 
culating the ratio TZ of the number of clusters with ^ 30° 
to the number with > 30°. This definition of the align- 
ment signal gives an indication of the steepness of the dis- 
tribution as well as the relative number of aligned and non- 
aligned cluster-BCG pairs. The errors are determined as- 
suming that the counts are Poisson processes with uncer- 
tainty a a = y/a where a is the count, and that the uncer- 
tainties in the two counts are independent. A random dis- 
tribution would yield TZ= 0.5. We show the distribution of 
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Figure 3. Comparison of the axis ratio b/a of BCGs and clusters. 
There is no correlation between the two axis ratios. 



Figure 5. Histogram of cluster-BCG alignment in the Total (solid 
line), Centre (dashed line) and Offset (dasli-dotted line) samples. 
The distributions are not fiat. The alignment signal for the Centre 
sample is approximately 5.5a stronger than in the Offset sample. 
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Figure 4. Distribution of member counts in clusters (Total sam- 
ple), which is used as a cluster-richness measure in the following 
analysis. The count includes red sequence galaxies up to 3 mag- 
nitudes fainter than the BCG. 

alignments in Figure |5] and tabulate TZ in Table |2] There is 
a clear and unambiguous alignment signal (the Binggeli ef- 
fect). In the Centre sample the alignment signal {TZ— 0.865) 
is approximately 5.5a stronger than in the OfTset sample 
{TZ= 0.659) where the significance is defined as the differ- 
ence between TZcentre and TZo f fset weighted by the uncer- 



tainties summed in quadrature. We find only a substantially 
weaker alignment for second ranked galaxies {TZ= 0.651 for 
centred second ranked galaxies, TZ= 0.533 for offset) and 
random alignment for third ranked galaxies. The observed 
ahgnment is unique to the BCG which is another character- 
istic setting it apart from other cluster galaxies. 

In the Offset clusters the line connecting the cluster cen- 
tre and the BCG is strongly aligned with the cluster shape, 
with TZ= 1.10 ± 0.04 (Figure [6|. In a bottom-up structure 
formation scenario with galaxies moving in a cluster that is 
oriented along a filament, a line drawn from a massive galaxy 
to the centre of the cluster will be strongly aligned with the 
direction of elongation of the cluster. However, we may in- 
stead be seeing a purely geometrical effect: in an elongated 
cluster more galaxies lie near the major axis of the cluster 
than the minor axis, hence a line drawn from most galaxies 
to the centre will be aligned with the cluster. To investi- 
gate this we also look at the alignment for the fainter red 
sequence galaxies farther than 200 kpc from the centre. If 
we are really seeing a purely geometric effect the alignment 
should be similar in these samples as well. We find that the 
samples have 7^= 0.94,0.92,0.92,0.86 for the second, third, 
fourth, and fifth ranked red sequence galaxies respectively. 
The alignment in the second ranked etc. galaxies is signifi- 
cantly non-random even if it is roughly 3a weaker than for 
the BCG. This suggests that the observed alignment is the 
result of a combination of influences. For the BCG, the ef- 
fect of movement along the filament is stronger than for the 
other galaxies, which is why TZ is higher. However, since all 
fainter galaxies show a similar alignment the geometry of 
the cluster is likely to have a significant influence as well. 
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Table 2. Alignment between BCG and host cluster. 



Dong et al. 



Sample 


^Cluster 
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Dong et al. 














(0.1 < z < 0.3) 






Kocster et al. 




Sample 


-^Cluster 


■R 


± 


-^Cluster 


n 


± 


Total 


2719 


0.833 


0.032 


5744 


0.858 


0.023 


Centre 


1365 


0.975 


0.052 


3647 


1.00 


0.033 


Offset 


1354 


0.710 


0.039 


2097 


0.649 


0.029 
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Figure 6. Histogram of alignment between the cluster and the 
BCG-centreline. The centreline is the line drawn from the cluster 
centre to the BCG. This test is only performed on the Offset 
sample whose BCGs are farther than 0.2 Mpc from the centre. 
We find a very strong alignment signal, with 1^= 1.10 it 0.04. 



3.3 The Effect of BCG Dominance 

One readily observed distinguishing property of BCGs is 
their dominance in luminosity over other cluster galaxies. 
They do not simply form the statistical extreme of the clus- 
ter galaxy luminosity function but require a distinct physical 
formation process from other cluster galaxies. A statistical 



test developed by Tremaine & Richstone ( 1977 1 and used by, 
among others, Loh & Strauss' (2006| has been used to show 



this using the difference between first and second ranked 
galaxies in clusters. The processes by which BCGs become 
luminous, however, are not well understood. Simulations by 
[Dubinski] (e.g. |1998[ ) predict that the BCG forms through 
mergers of several massive galaxies moving along a filament 
early in the cluster's history. Since the alignment effect may 
also result from a preferred infall direction of galaxies dur- 
ing cluster formation a connection between dominance and 



alignment is worth investigating. Figure [7] shows the distri- 
bution of the BCG's dominance, in the Dong et al. sample, 



defined following Kim et al. ( 2002 1 as the difference in 



band magnitude of the BCG and the mean magnitude of 
the second and third ranked galaxy. 



dom ; 



; nil 



m2 + ma 



(9) 



An alternative would be to investigate the dependence 
of alignment on the absolute magnitude of the BCG. How- 
ever, if the BCG grows by merging with the second and third 
ranked galaxies this measure of dominance is more appro- 
priate. The dominance distribution is more strongly skewed 



than that of Kim et al. ( 2002 \ , presumably because of the 



differences in the details of the cluster selection algorithms. 
The mean dominance for the Centre sample is 0.90 mag and 
the mean dominance for the Offset sample is 0.81 mag. |Lohl 
fc Strauss] ( [2006 [ ) find a magnitude difference between first 
and second ranked galaxies of 0.8 mag, which is in good 
agreement with our values. 

If the BCG's high luminosity and its alignment with 
their host cluster stem from infall and merging along a pre- 
ferred direction one might expect a correlation between BCG 
dominance and its degree of alignment with the cluster. [Km] 
et al.j ('2002 ) found in a small SDSS sample that 39 out of 66 
dominant BCGs show strong alignment (A^ < 30°) but find 
no statistically significant alignment in a sample of 49 non- 
dominant BCGs. In a scatter plot of dominance vs. align- 
ment (Figure [8^) we already see an indication of stronger 
alignment (smaller A(/!>) in high dominance clusters. Defin- 
ing as dominant those BCGs that have domjs 0.65, which 
corresponds roughly to the peak of the dominance distribu- 
tion (Figure [7| we find that the dominant sample has an ap- 
proximately 4.4(T stronger alignment than the non-dominant 
sample (see Figure [sJd and Table ]3] for TZ values). However, 
the non-dominant sample is not completely isotropic. In the 
Centre cut the dominant sample is more strongly aligned 
than the non-dominant sample (approximately 3.7a). The 
difference between the dominant and non-dominant sub- 
sample of the Offset cut is significant at only the 2.1a" level. 

Using Kolmogorov-Smirnov tests we find with high con- 
fidence (Table ]4| that the A(f) distribution of the total dom- 
inant sample is drawn from a different underlying distribu- 
tion than the non-dominant sample. Similarly we find that 
the Centre dominant sample comes from a different underly- 
ing distribution than the non-dominant sample. In the Off- 
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Figure 7. Histogram of BCG dominance over the other cluster 
galaxies. Dominance is defined as the difference between the i- 
band magnitude of the BCG and the average of the second and 
third ranked galaxies. In the analysis of the alignment effect, we 
split the sample into dominant (dom> 0.65) and non-dominant 
BCGs. 



Table 4. Kolmogorov-Smirnov tests determining whether the 
dominant and non-dominant sub-samples are drawn from differ- 
ent underlying distributions. A gives the significance of the dif- 
ference between TZ statistics for the dominant and non-dominant 
samples in units of standard deviation (cr). 



Sample 


K-S Confidence 


A 


Dong Total (0.08 < z < 0.44) 


99.99% 


4.4 


Centre 


99.92% 


3.7 


Offset 


96.52% 


2.f 


Dong Total (O.f < z < 0.3) 


99.99% 


3.4 


Centre 


99.98% 


3.0 


Offset 


82.92% 


1.3 


Koester Total 


99.99% 


4.7 


Centre 


99.76% 


3.1 


Offset 


9f.39% 


2.4 



set sample the difference between the dominant and non- 
dominant samples is less significant. 

3.4 The Effect of Cluster Richness 

The denser environments of rich clusters offer more possibil- 
ities for galaxy-galaxy interactions during cluster formation 
and cluster evolution,which may affect the observed align- 
ment. Figure|9^ is a scatter plot of richness i.e. number of red 
sequence galaxies in a given cluster with alignment. There is 
a slight indication that richer clusters show stronger align- 
ment. Figure [9|3 shows the histogram of alignment for clus- 
ters with fewer than and more than 20 members, picking 
out the high-count tail of the cluster richness distribution 



Table 6. Kolmogorov-Smirnov tests determining whether the rich 
and poor sub-samples are drawn from different underlying dis- 
tributions. A gives the significance of the difference between TZ 
statistics for the rich and poor samples in units of standard devi- 
ation (cr). 



Sample 


K-S Confidence 


A 










qq qq% 


2.3 


Centre 


99.88% 


2.2 


Offset 


92.04% 


0.6 


Dong Total 






(0.1 <z < 0.3) 


97.94% 


1.4 


Centre 


90.89% 


1.2 


Offset 


77.37% 


0.7 


Koester Total 


96.53% 


1.7 


Centre 


96.88% 


1.6 


Offset 


43.98% 


0.7 



(Figure |4|. Rich clusters show a stronger alignment, a dif- 
ference significant at the 2.3(t level. The rich clusters in the 
Centre subsample are even more strongly aligned. With the 
offset cut, the results are not significant (Table [5|. K-S tests 
show that the rich and poor samples in the Total and Centre 
cuts are different with high confidence, but the difference is 
statistically insignificant in the Offset sample (see Table [6|. 

We investigate the possibility that position angle de- 
terminations may be noisier in poor clusters than in rich 
clusters due to small number statistics which may suppress 
the alignment signal. This could explain the weaker align- 
ment seen in poorer clusters. In Figure [lO] we show the sta- 
tistical uncertainty on the cluster position angle (cr^) as a 
function of cluster richness (left-hand panel) and cluster axis 
ratio (right-hand panel). The main source of uncertainty in 
our determined cluster position angles stems from a clus- 
ter's axis ratio. With this in mind we investigate the axis 
ratio distributions of our rich and poor sub-samples. The 
mean axis ratio in the rich sub-sample is higher (i.e. the 
clusters are rounder) {b/a) = 0.53 than that in the poor 
sub-sample {b/a) — 0.38. Since a higher axis ratio will de- 
grade the alignment signal more by introducing larger shot 
noise errors in the cluster's position angle, we are confident 
that the stronger alignment observed in the rich clusters is 
real. 

Even though rich clusters tend to have more strongly 
aligned BCGs than poor clusters, their BCGs are less dom- 



inant than those of the poor clusters (Figure 111. Assuming 
no correlation between cluster richness and the galaxy lumi- 
nosity function (see e.g.. [Tremaine fc Richstone|[T977[ |Loh| 
fc Strauss||2006 1 it is more likely that a rich cluster has an 
uncommonly bright second ranked cluster galaxy, resulting 
in less dominant BCGs than in poor clusters. 



3.5 Redshift Evolution 

It is certainly possible that the mechanisms responsible for 
the alignment between BCGs and the cluster act over the 
full span of the cluster's lifetime rather than during or imme- 
diately after cluster virialization. For example, tidal torques 
exerted by the cluster or secondary infall may enhance or 
reduce any primordial alignment (e.g.. Ciotti fc Dutta|1994 



© 2010 RAS, MNRAS 000,[T]{T5] 



Alignment of Brightest Cluster Galaxies with their Host Clusters 9 



0.20 



0.15 
0.10 
0.05 



0.00 



100 




3.0 



(a) 



20 40 60 80 100 
A<f> (Deg.), Dominont Somple 



0.20 



0.15 = - 



0.10 



0.05 



0.00 




20 40 60 80 100 
A0 (Deg.), Not Dominont Somple 

(b) 



Figure 8. Panel (a): Scatter plot of dominance vs. alignment. There is a clear indication that dominant clusters tend to be more aligned 
than non-dominant clusters. Panel (b): Histogram of cluster-BCG alignment examining the effect of BCG dominance on the alignment 
signal. The top panel shows the dominant sample (dom^ 0.65), and the bottom panel shows the non-dominant sample. The Total (solid 
line) and Centre (dashed line) samples have a significant difference between dominant and non-dominant clusters. The Offset sample 
shows a smaller, statistically less significant, difference between the two sub-samples. The statistics for these various samples are given 
in Table |3] 



Hopkins et al. 



Merritt ( 1985 1 and 



2005 



Altay et al.|2006 1 . It has been argued by 



Tremaine ( 1990 1 among others, that the 



observed dominance of BCGs cannot be achieved via can- 
nibalism of other cluster members during the BCG's life- 
time, since the high velocity dispersion of cluster galaxies 
makes frequent merging unlikely. Thus the BCG's domi- 
nance must already be in place shortly after the cluster has 
formed. Indeed, we observe no evolution in the distribution 
of BCG dominance from high to low redshift (at least for 
0.08 < z < 0.44). Thus if we observe any evolution of align- 
ment during the cluster's lifetime it cannot be the result 
of the same mechanism that makes the BCGs dominant. 
To test this, we split our dominant and non-dominant sam- 
ples into a high redshift and a low redshift sample, taking 
z — 0.26 roughly the median of the Dong et al. sample as the 
dividing line. The median look-back times of the two sub- 
samples are 2.1 Gyr and 3.7 Gyr. We find that the mean axis 
ratios of the BCGs and the clusters are the same in the high 



and low redshift samples. Thus uncertainties in measuring 
the BCG's shape stemming from lower signal-to-noise ratio 
photometry and a smaller physical size on the sky at high 
redshift, do not bias our alignment signal. 

Figure [12] shows the alignment histograms for the Total 
dominant and non-dominant samples at high and low red- 
shift. The alignment effect in the Total sample increases in 
TZ by approximately 3. la from high to low redshift. In the 
dominant sample we see an increase of roughly 2.4cr; there 
is no statistically significant evolution in the non-dominant 
sample. Even though the TZ values for the high and low red- 
shift samples are only 2a apart, K-S tests show at the 99% 
confidence level that the high and low redshift samples are 
drawn from different underlying distributions in both the 
dominant and non-dominant samples for the Centre sample. 
The TZ statistic probes a specific aspect of the alignment 
distribution where as the K-S tests asks the more general 
question of similarity between distributions. Two different 
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Table 3. Alignment between BCG and clusters for dominant and non-dominant BCGs. 



Dong et al. 
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^Cluster 


n 


± 








Total Dominant 


4353 


0.818 


0.025 








Total Not-Dominant 
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n 


± 
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7^ 
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Total Dominant 


1692 


0.921 


0.045 


4009 


0.927 


0.029 


Total Not-Dominant 


1025 


0.705 


0.045 


1735 


0.715 


0.035 


Centre Dominant 


907 


1.090 


0.072 


2722 


1.065 


0.041 


Centre Not Dominant 


457 


0.778 


0.073 


924 


0.844 


0.058 


Offset Dominant 
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0.756 


0.055 
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0.039 


Offset Not Dominant 
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Figure 9. Panel (a): Scatter plot of richness vs. alignment. Rich clusters show a stronger preference for alignment than poor clusters. 
Panel (b); Histogram of cluster-BCG alignment examining the effect of cluster richness on the alignment signal. The top panel shows 
the rich samples (Total sample solid line, Centre sample dashed line, Offset sample dash-dotted line) (A'^Members ^ 20), and the bottom 
panel shows the poor samples. The alignment signal in the total rich sample is approximately 2.3(t stronger than in the poor sample 
(Table ig [el. 
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Figure 10. Left panel: Position angle uncertainty as a function of cluster axis ratio. The position angle uncertainty is determined by 
propagating the shot- noise errors on the Stokes parameters through Equation |6] Right panel: Position angle uncertainty as a function of 
cluster richness. There is a clear trend towards higher position angle uncertainty in rounder clusters. The effect of the number of cluster 
members is much less clear. 




(a) (b) 

Figure 11. Panel (a) Scatter plot of cluster richness and BCG dominance, indicating that poorer clusters have more dominant BCGs. 
Panel (b) Histogram of BCG dominance split on our definition of rich and poor clusters. The BCGs of rich (solid line) clusters tend to 
be less dominant than those of poor (dashed line) clusters. In a rich cluster it is more probable that there is an unusually bright second 
ranked cluster galaxy than in poor clusters, assuming that both cluster galaxy populations have the same luminosity function. 
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Table 5. Alignment between the BCG and cluster for rich {N Member ^ 20) and poor {N Member < 20) clusters separately. 



Dong et al. 



Sample 


^Cluster 


n 


± 








Total Rich 


259 


1.023 


0.120 








Total Poor 


6772 


0.744 


0.018 








Centre Rich 


141 


1.350 


0.229 








Centre Poor 


3272 


0.849 


0.030 








Offset Rich 


118 


0.735 


0.137 








Offset Poor 


3500 


0.656 


0.023 
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Ivoester et al. 




Sample 


^Cluster 


n 


± 


^Cluster 


■R 


± 


Total Rich 


81 


1.190 


0.265 


211 


1.110 


0.153 


Total Poor 


2638 


0.823 


0.032 


5533 


0.849 


0.023 


Centre Rich 


41 


1.563 


0.500 


135 


1.368 


0.238 


Centre Poor 


1324 


0.961 


0.053 


3512 


0.992 


0.033 


Offset Rich 


40 


0.905 


0.286 


76 


0.767 


0.178 


Offset Poor 


1314 


0.704 


0.039 


2021 


0.644 


0.029 



distributions can have the same Ti, value but a K-S test will 
show that they are difTerent. However, due to Poisson noise 
two similar distributions may have different TZ values yet the 
K-S test will show that they are in fact the same. Thus in 
this case it seems that even though the TZ values are only 
2(T apart the results of the K-S test are a stronger constraint 
on the redshift evolution. The statistics are summarised in 
Tables [7|8l andflOl 

If we now split the Total sample on richness instead of 
dominance (Tables |9] and 10 1, we find no statistically sig- 
nificant evolution in the alignment of the rich sub-sample. 
The Poisson errors are substantial, however, due to the small 
number of rich clusters in the sample. However, there is a 
trend toward stronger alignment at low redshift in the poor 
sample, which is significant at 3.2(j. 



4 SUMMARY AND DISCUSSION 

Using the SDSS Data Release 6 and the cluster catalogues 
develo ped by [Dong et al. ( 2008) (7031 clusters) and |Koester| 
|et al.| ( |2007b| (5744 clusters) we investigate the alignment 
between the BCG and cluster position angles, searching 
for infiuences from BCG dominance in luminosity, cluster 
richness, and redshift. Understanding the alignment offers 
an observational test of the current cosmological paradigm 
and also allows for a better understanding of contamination 
in lensing studies ( [Mandelbaum et al.| |2006[ ) . We have 
confirmed that BCGs are preferentially aligned with their 
host cluster and find that dominant BCGs have a AAa 
stronger alignment signal than non-dominant ones. This 
suggests that dominance and alignment result from the 
same physical process. The sample of clusters in which 
our selected BCG is farther than 0.2 Mpc from the cluster 
centre shows a 5.5cr weaker alignment signal, which may 
be an indication that these Offset clusters are the result 
of cluster mergers. If dynamical interactions influence 
alignment, one might expect a dependence of alignment 
on cluster richness. We found that rich clusters (defined as 
containing more than 20 red sequence galaxies within 0.5 
Mpc from the cluster centre) show a stronger alignment 



Table 10. Kolmogorov-Smirnov tests determining whether the 
high-redshift and low-redshift sub-samples are drawn from differ- 
ent underlying distributions. A gives the significance of the dif- 
ference between "R. statistics for the high and low redshift samples 
in units of standard deviation (cr). 



Sample 


K-S Confidence 


A 


Total 


99.14% 


3.1 


Centre 


99.82% 


3.2 


Offset 


38.25% 


1.0 


Total Dominant 


97.75% 


2.4 


Total Not Dominant 


80.45% 


1.7 


Total Rich 


11.91% 


0.4 


Total Poor 


99.48% 


3.2 


Centre Dominant 


99.07% 


2.1 


Centre Not Dominant 


98.97% 


2.4 


Centre Rich 


36.80% 


0.3 


Centre Poor 


99.83% 


3.3 


Offset Dominant 


46.37% 


1.1 


Offset Not Dominant 


49.27% 


0.1 


Offset Rich 


36.08% 


0.2 


Offset Poor 


38.66% 


1.0 



signal than do poor clusters. The difference in alignment 
is only significant at the 2.3cr level, but this is due to the 
large Poisson errors given the small sample of rich clusters. 
A Kolmogorov-Smirnov test shows that the alignment 
distributions in the rich and poor clusters are different at 
better than 99% confidence. 

Our cluster sample extends out to 2 = 0.44, allowing us 
to test redshift evolution. Non-linear effects during cluster 
virialization can weaken primordial alignments. Thus if the 
BCG has resided in the cluster for many crossing times its 
initial infall direction may have been forgotten. However, 
we find that the alignment signal is stronger at lower 
redshift, an effect significant at more than 3(t for those 
clusters whose BCGs lie within 0.2 Mpc of the Centre. For 
the Offset clusters we do not find a statistically significant 
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Table 7. Alignment between the BCG and cluster distinguishing between high redshift (z ^ 0.26) and low redshift {z < 0.26) clusters. 







z ^ 0.26 






z < 0.26 




Sample 


^Cluster 


n 


± 


^Cluster 


n 


± 


Total 


5036 


0.718 


0.021 


1991 


0.851 


0.038 


Centre 


2418 


0.803 


0.033 


992 


1.037 


0.066 


Offset 


2618 


0.646 


0.026 


999 


0.696 


0.045 



Table 8. Alignment between the BCG and cluster distinguishing between dominant and non-dominant BCGs at high and low redshift. 







z > 0.26 






2 < 0.26 




Sample 


^Cluster 


n 


± 


^Cluster 


n 


± 


Total Dominant 


3077 


0.779 


0.028 


1276 


0.919 


0.051 


Total Not Dominant 


1959 


0.630 


0.029 


715 


0.740 


0.056 


Centre Dominant 


1567 


0.890 


0.045 


689 


1.082 


0.082 


Centre Not Dominant 


851 


0.662 


0.046 


303 


0.942 


0.108 


Offset Dominant 


1510 


0.676 


0.035 


587 


0.757 


0.063 


Offset Not Dominant 


1108 


0.606 


0.038 


412 


0.616 


0.062 



evolution in alignment over the redshift range considered 
here. 

A likely cause for the correlation between BCG 
dominance and its alignment with the host cluster galaxy 
distribution is that both are the result of galaxies falling 
along primordial filaments and merging during cluster 
formation (see e.g. 



West et al. 



1995 



Fuller et al. 



1999 



Torlina et aL] |2007| [Garijo et al.||1997| |Dubinski| [T998 
Knebe et al.|2004 i. Recent modelling using the Millennium 
Simulation ( |De Lucia fc Blaizot [2007 1 suggests that the 
stars that eventually end up in the BCGs are formed early 
{z « 4) but that the final BCG is assembled comparatively, 
late with half their mass being locked up in a single galaxy 
only after z « 0.5. The stellar mass of BCGs grows by a 
factor of 3 — 4 via mergers since z = 1 (Gao et al. 2004 



Aragon-Salamanca et al. 1998 De Lucia & Blaizot 20071 



Observational studies, however, suggest that 80% of the 
stars that make up massive (> 4L*) red galaxies today are 
bound up in a single galaxy hy z = 0.7 and that the stellar 
mass in red galaxies evolves only by a factor of 2 since 
2 = 1. In the case of massive red galaxies this is mostly 
due to dry mergers ( |Faber et aL]|2007| [Brown et al.||2007 |. 
The observation that BCG dominance evolves little with 
redshift (at least in the range considered here) suggests that 
the mechanism that causes the BCG to become dominant 
acts early on in the cluster's history. However, as pointed 



out by Faltenbacher et al. (20071 any primordial alignment 



on small scales is likely reduced by non-linear processes 
during cluster virialization. A possibility is that the BCG 
is somehow resistant to realignment after it has settled at 
the bottom of the potential well. Secondary infall episodes 
may reinforce the primordial alignments, which may be the 
cause of the slight redshift evolution of the alignment that 
we observe. If galaxy mergers early in a cluster's history 
are the cause of both BCG dominance (with more mergers 
creating more dominant BCGs) and its alignment with the 
host cluster, the uniformity of the BCG population as a 
whole remains puzzling, since BCGs in poorer clusters will 
have undergone fewer mergers than those in rich clusters. 
Perhaps the colour and luminosity of BCGs are not affected 
by mergers after a certain threshold of merging activity. 



Additional mergers then only affect the BCG orientation 
within the cluster removing other galaxies from the cluster's 
population (preferentially brighter more massive galaxies) 
letting the BCG grow in dominance. Further investigations 
of clusters in the redshift range 0.5 < z < 1 will be very 
useful in understanding how the growth in dominance of 
the BCG and its alignment are related. 

Filamentary infall is likely to play an important role in 
determining BCG and cluster orientations early in cluster's 
history, but tidal torques from the large scale cluster mass 
distribution can influence the alignment signals between 
BCGs and clusters during the entire cluster lifetime. Sim- 
ulations show that tides align a galaxy's major axis with 
the cluster's radial direction on timescales greater than 
the galaxy's dynamical time scale but significantly shorter 
than a Hubble time. The effect is strongest for galaxies 
outside the cluster core radius ( Ciotti fc Dutta|1994 1 , and is 
therefore less important for central BCGs. The simulations 



by Faltenbacher et al. (20081 find on galaxy scales tides 



cause a radial alignment. Even though tidal torques cause 
some kinds of alignment they are not likely the source 
of the BCG-cluster alignment we observe. The evolution 
of the alignment signal needs to be better understood 
theoretically to differentiate the effects of tidal torques and 
infall episodes. These questions can likely be addressed 
using the Millennium Simulation which has already been 



used to explore properties of BCGs (see e.g. Ruszkowski & 
Springel|[2509l ) 



Recent surveys such as the SDSS Stripe 82 and future 



surveys such as LSST (LSST Science Collaborations: Paul 
|A. AbeU et al.l2009| LSST Science Book) and PanSTARRS 
( Kaiser et al.| 2002 \ with deeper photometry will allow 



us to probe clusters at higher redshifts to study the 
evolution of alignment over a larger range of cosmic time. 
Deeper surveys will also increase the number of galaxies in 
low-redshift clusters, reducing Poisson noise in our samples. 
This will allow us to study cluster shapes and alignments 
for galaxies of different luminosities, a diagnostic of different 
accretion and dynamical histories. Simulations addressing 
the connection between BCG dominance and alignment will 
be important in discovering which mechanisms cause the 



© 2010 RAS, MNRAS OOP, [1^151 



14 M.Niederste-Ostholt et al. 



Table 9. Alignment between the BCG and cluster distinguishing between rich (Af^v/embei s ^ 20) and poor (Nj^i^^i,^^ < 20) clusters at 
high and low redshift. 







z ^ 0.26 






z < 0.26 




Sample 


^Cluster 


7^ 


± 


-^Cluster 


■R 


± 


Total Rich 


203 


1.051 


0.148 


56 


0.931 


0.249 


Total Poor 


4836 


0.706 


0.021 


1936 


0.847 


0.039 


Centre Rich 


112 


1.383 


0.265 


29 


1.231 


0.460 


Centre Poor 


2308 


0.782 


0.033 


964 


1.029 


0.066 


Offset Rich 


91 


0.750 


0.159 


27 


0.688 


0.269 


Offset Poor 


2528 


0.642 


0.026 


972 


0.696 


0.045 
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A0 (Deg.), Domlnont Somple 




20 40 60 80 100 
(Deg.), Not Dominant Sample 



Figure 12. Histogram of cluster-BCG alignment examining the 
effect of BCG dominance and redshift on the alignment signal. 
The top panel shows the dominant sample and the bottom panel 
shows the non-dominant sample. The solid histograms show the 
high redshift (z ^ 0.26) sample and the dashed histograms show 
the low redshift sample. The alignment increases by approxi- 
mately 2.4(7 from high to low z in the dominant case. In the 
non-dominant case there is no significant change in alignment. 



distinct physical properties of the BCGs and their relation 
to their parent cluster. 
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